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We analyze the possibility of detection of cosmological relic neutrinos via neutrino capture on beta 
decaying nuclei. This reaction has no threshold in neutrino energy, which is crucial for searching for 
relic neutrinos possessing very low energies. We focus on tritium (^H) and rhenium (^*^Re) beta 
radioactive isotopes to be used in KATRIN and MARE experiments dedicated to measurement of 
the electron neutrino mass at sub-eV scale. We examine these experiments from the viewpoint of 
searching for the cosmological neutrinos via neutrino capture. We conclude that even with possible 
gravitational clustering of relic neutrinos the prospects for their detection in these and other similar 
experiments are not optimistic. Nevertheless KATRIN and MARE experiments could establish some 
usefull constraints on density of relic neutrinos. 
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I. INTRODUCTION 

The cosmological relic neutrinos are one of the most important and abundant constituents of the Universe predicted 
by Big Bang cosmology. There are about 10®'' neutrinos per flavor in the visible Universe [l|, which corresponds to the 
number density per flavor of relic (anti-) neutrinos in average over the Universe about (ry) ~ 56 cm~^0. In number, 
neutrinos exceed the constituents of ordinary matter (electrons, protons, neutrons) by a factor of ten billion. On the 
other hand their existence has not yet been confirmed by direct observations. This represents a challenging problem 
of modern cosmology and experimental astroparticle physics. Many proposals in the literature aim to observation of 
indirect astrophysical manifestations of relic neutrino sea However laboratory experiments would be the most 

robust probe of this component of the Universe. There are various strategies discussed for laboratory searches of relic 
neutrinos [10-20.]. The neutrino capture reaction, 19] 

^e + {A,Z)^{A,Z + l) + e-, (1) 

is one of the most discussed possibilities in the literature 0, [l7l - l20l | . Since relic neutrinos have very low energy this 
reaction should have no threshold implying that the initial nucleus is /3-radioactive. This idea was first proposed many 
years ago by Weinberg in Ref. @ , but in the context of massless neutrinos with large chemical potential. The latter 
proved to be inconsistent with the Big Bang cosmology. The discovery of neutrino mass revived the interest to 
the above reaction, since at the endpoint of the electron spectrum the energy difference between /3-decay and neutrino 
capture electrons should be around ~ 2m^. In principle, this effect can be observed with a detector having an energy 
resolution less than this difference. However, the existing experimental techniques are unable to detect relic neutrinos 
via neutrino capture if their number density in our vicinity does not significantly differ from the global average number 
density (rj) ^ 56 cm~^ [l8l - l20| . As it was pointed out recently ilSi, i21, ^22J , the gravitational clustering of neutrinos in 
our galaxy or galaxy cluster may drastically enhance the relic neutrino density in comparison to the overall average 
{rj) making its detection more realistic. 

In our previous papers [l^ H^l we studied reactions of neutrino capture on single ([T]) and double beta decaying 
nuclei as a tool for probing cosmological relic neutrinos. Here we focus on detailed calculation of the event rate of the 
reaction ([T]) and estimate typical detection rates of the relic neutrino capture taking into account possible effect of 
gravitational clustering. In Sec. |TT]we present our approach to nuclear /3-decay with an emphasis on tritium '^H and 
rhenium ^^''Re radioactive isotopes. In Sec. Illll we apply these approach to analysis of relic neutrino capture on these 



isotopes and discuss KATRIN '^H [23, [2J] and MARE ^*^Re experiments 25| dedicated to measurement of electron 
neutrino mass at sub-eV scale. We examine prospects of these experiments for searching for the relic neutrinos via 
neutrino capture on the corresponding nuclei. 

II. NUCLEAR /3-DECAY. BASIC INGREDIENTS 

Here we discuss theoretical aspects of single /3-decay of nuclei with special emphasis on tritium and rhenium isotopes 
[26j . We introduce nuclear structure ingredients used in the subsequent sections for analysis of relic neutrino capture 
by these nuclei. 

Let us consider general case of nuclear /3-decay 

{A,Z)^ {A,Z + l)+e- (2) 
driven by the standard weak /3-decay Hamiltonian 

n^{x) = ^e(x)7'^(l - jMx)j^{x) + h.c. (3) 

Here, = GpcosOc, where cos^c is the Cabbibo angle. e{x) and i^e{x) are the electron and neutrino fields, 
respectively. The strangeness conserving free nucleon charged current is 

jf'ix) = p(x)7^(5y - gA75)n{x), (4) 

where p{x) and n{x) are the proton and neutron fields, respectively, gv — 1.0 and gA — 1.25. 

Single /3-decay occurs in the first order in the weak interaction. The corresponding S-matrix element is given by: 

= 2Ti5{Ef +Ee + E^- E,)i-t)—^ X 

v2 

Mx, E,)-ff'il - j5WAx, E,){A, Z + 1| J^(0, x)\A, Z)dx. (5) 
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Here, is the nuclear weak charged current in the Heisenberg representation. In our notation Ei, Ef, E^ and E^, are 
energies of initial and final nuclei, electron and antineutrino, respectively. The wave functions of outgoing electron 
and antineutrino are denoted as ipei'^Pt- We use non-relativistic impulse approximation for the nuclear current: 

A 

J^(0, x)^Yl ^« [5y5^° + 9A{<yk)ng^'']S{x - x^). (6) 

n=l 

From now on it is convenient to consider /3-decay of tritium and rhenium separately. 



A. /3-decay of tritium 



In the tritium /3-decay the spin and parity of initial and final nuclei are equal: 

3iJ((l/2)+) ^ 3^e((l/2)+) + e- + F„ 



(7) 



In this case the dominant contribution to the decay rate is determined by Si/2 wave-states of outgoing electron and 
antineutrino 



Mx,Ee) « ^Fo{Z + l,Ee) u{P,) (8) 
rAx,E,)^u{~P,,) (9) 

where x is the coordinate of the lepton. Fq is the relativistic Fermi function (for more details see Appendix) J27']. 
Normalization of the spinor is u^{P)u{P) = 1. Pe = {Ee,Pe) and = [E^^p^) are 4-momenta of electron and 
antineutrino, respectively. 

Substituting Eqs. (jS]), © and © into Eq. ([S]) we get for T-matrix clement an expression 

(/|rW|z) = H)-^V^o(^ + i,i^e) X 

where Fermi and Gamow- Teller nuclear matrix elements are defined as 

Mj.(to,to') =3ae((l/2)+m'|^7:+|(l/2)+TO)3H, 



(10) 



MGT(m,TO')=^//e(l/2+™'|^r+<T„|(l/2)+m)3^. 



(11) 



Here, m,m' = ±1/2. 

Now we sum up over the lepton polarizations, over the projection of angular momenta, m', of final nucleus and 
average over projection of angular momenta, m, of initial nucleus. Evaluating the corresponding traces for the squared 
T-matrix we have 

^ \{f\T('^^)\' = 2 (^^y Fo{Z + l,E,){BFCH)+BGTeH)) (12) 



with 



m,m' 



(13) 



Tn.m' 



We neglected terms proportional to the space components of the lepton momenta as they vanish after integration over 
angles in calculation of the total /3-decay rate. 

The Fermi and Gamow- Teller beta strengths in Eq. (jl3p can be represented in terms of reduced nuclear matrix 
elements: 

2 



>He((l/2)+ llE^n II (V2) + )3^ 



.^e((l/2)+ ||^r,>„|l (l/2)+)3^ 



(14) 
(15) 
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With our normalization of the lepton wave functions ([8]), ^ the differential decay rate is 

dT^ = Y: \{f\T^'W ^-S{E. + Ee + Ef - E.) ^ (16) 

Inserting Eq. (|12l) for squared T-matrix element and integrating over the phase space we find the total /3-decay rate 
for the reaction ([7]) 

^^{'H) = ^m, {Gpmlf {Bf{^H) + Bgt{'H)) I^i^H) (17) 
with the phase space integral 

I^i^H) = ^ / " ' Fo{Z + 1, Ee)peEe{E, - Ef - E.fdE,. (18) 

The dependence of H) on small neutrino mass is ignored. For tritium numerical integration in (|18p gives 

/'^(^i/) = 2.88 X 10"'^. (19) 

The nuclear matrix elements for tritium are known in the literature. We take for these matrix elements the values 
|MfP = 1, |MgtP = 3 derived in Ref. [28,]. Then we find 



which is very close to the measured value of tritium /3-decay half-life. Vice versa, using Eqs. ((T7| - p9|) . one can extract 

1/2 



a value of |A/i?|^ + g\ \Mgt\^ in a model independent way from the experimental value T^.^i^H) « 12.33 y. In this 



way we obtain 

IMflVgilMcTl' « 5.645. (21) 
This model independent estimate we will use in discussion of relic neutrino capture by tritium. 

B. /3-decay of rhenium 

Single /3-decay of rhenium 

i«^i?e((5/2)+) ^ i»^0s((l/2)-) + e- + V,. (22) 

is more complicated than the tritium /3-decay mainly because of the fact that the initial and final nuclear spins and 
parities are different. In (29| it was shown that the dominant contribution to this decay is given by the emission of 
electron and antineutrino in and Si/2 wave states, respectively. The wave function of electron has the form 

^P^f^{x,E,)c^i ^Fi{Z + l,E,) (^a;-pe + ^7-a=7-Pe^ u{P,). (23) 

Fi is the relativistic Fermi function (for more details see Appendix) (27| . 

Now, we are to analyze the squared matrix element in Eq. ([5]) with the electron and neutrino wave functions given 
in Eqs. (j!]), (|23|) . Towards this end we consider the quantity 

M = ^-)e{x, Ee)-/f,{l - -/5)'ip^{x, Ey) X 

((l/2)-m'|^r+(gv5^0-F5^g^'^(^„)'=)|(5/2)+™) 5{x-x^), (24) 

n 

which is the integrand of the space integral in Eq. ([5]). Carrying out summation over polarizations of leptons we keep 
only the terms surviving after integration over phase space. Then, for the squared T-matrix element we obtain [26j 

\{f\T^%)\^ = Gl g\ Mg:)'^^Fi(Z + l,i?,). (25) 



5 



where the squared nuclear matrix element includes summation over all the spin orientations of the final nucleus m! 
and averaging over the spin orientations of the initial nucleus m. The /3-strength for ^^'^ Re takes the form 



BGTrRe)=g\ 



M. 



(2-) 
GT 



1 



6 i?2 



((l/2)-||^r+{<T„$5r„}2 II (5/2)+) 



(26) 



Inserting the above expression to Eq. (|16p and integrating over the phase space of final electron and neutrino we 
find the total /3-decay rate of rhenium 



r^(i«7i?e) = "le^^g^ ^^^^ BcA^^'Re) I^C^'Re), 



where the phase space integral is 



I^{^»^Re)^^ /" ' ' p,E, {E,^Ef-E,fFi{Z + l,E,)pldE,,. 

We neglected the effect of small neutrino mass. For rhenium numerical integration in ([28]) gives 

/'3(i87^g) ^ 1 22 X 10-^. 

With this value of the phase space integral and using the experimental half-life 

rf/2(^*^^i?e) = 4.35 X 10^° years 
we obtain the value of rhenium nuclear matrix element 
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M. 



(2-) 
GT 



3.57 X 10" 



(27) 



(28) 



(29) 



(30) 



(31) 



We will use this value in our analysis of neutrino capture on rhenium where the same quantity will appear. 



III. NEUTRINO CAPTURE WITH BETA DECAYING NUCLEI 



Now let us consider processes, similar to the previously considered /3-decays, but with neutrino in the initial state 
of reaction: 



(32) 



This neutrino can be of any origin. An interesting possibility is a neutrino from the cosmic neutrino background. 
Thus the above reaction may represent one of the very few possible ways to detect this important component of the 
Universe. Since the energy of cosmic neutrinos is negligibly small, the above neutrino capture reaction takes place 
only if there is no threshold. This implies that the target nucleus should be /3-radioactive. In the previous sections we 
considered in details two candidates of this type, tritium and rhenium. Below we study capture of cosmic neutrinos 
with these radioactive nuclei. 

The differential rate of the reaction (15^ is given by 



dpe 

(27r) 



(33) 



for the normalization of incoming neutrino plane wave to be one particle per volume V. 

The kinetic energy of electron produced in /3-decay is continuously distributed within the interval < E^ — rrie < 
Qp ~ rrii, while in the reaction ([5^ the final electron is monoenergetic with E^ ~ rrie = Qp + El,'^\ where £'^^' is the 
relic neutrino energy. Thus the signatures of /3-decay and cosmic neutrino capture are quite different. However the 
energy difference between /3-decay and neutrino capture emitted electrons is very small ^ 2mi, which is challenging 
for energy resolution of detectors used in experiments searching for this process. ^From the viewpoint of matrix 
element calculations neutrino capture reaction ((32)) with very low energy neutrinos, such as the cosmic neutrinos, is 
very similar to /3-decay of the target nucleus. 

Below we consider two target nuclei tritium and rhenium which are planned to be used in the experimental setups 
KATRIN with ■^H and MARE with ^^^Re aimed to measure neutrino mass in sub-eV range. 
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A. Neutrino capture on tritium 

The cosmic neutrino capture reaction by tritium, 

V, +3 i/((l/2)+) ->3 i/e((l/2)+) + e-, (34) 

is characterized by the T-matrix element of the same form as in the case of single /3-decay. Within the same approx- 
imations and summing over the polarizations we have 

^ |(/|T(i)|*)|2=G^Fo(Z + l,i?e)(i?F('i?) + SGT('i/)). (35) 
Inserting (f35|) to ([33l) . for the total capture rate we find 

V''(^H)^\--GlF^{Z + \,E,){BFi^H)+BGTQll)) E,, (36) 

V TT 

where = Ei + Ei, — Ef. With our neutrino wave function normalization both expressions correspond to capture 
rate per neutrino per '^H atom. For the number density of cosmic neutrinos rj^ we replace 1/V — > rji, in (j36p and get 
the capture rate per atom 

T-{^H) - -Gl Fo{Z + 1,E,) E, (BfCH) + BgtCH)) (r;,). (37) 

Here rj^, is the local cosmic neutrino number density which could be significantly larger than the average over the 
universe (rji,) ^ 56 cm~'^ due to gravitational clustering [Tsl. [2TI| 

-l3 1 n4 



W - 10^ (38) 



assuming tUi, — leV and baryon density rjt, = 10^^ — lO^^'cm^'^ in a cluster of galaxies. 

A combination of the nuclear matrix elements in Eq. (1371) . shown in curl brackets, we previously encountered in 
/?-decay rate and found its numerical value from the half-life T1/2 of "^H. The corresponding value is given in Eq. (j2ip 
which we substitute in ([37]) and obtain 

r''(3i/) =4.2 X 10-25^ y-i. (39) 

This value is compatible with the result of Ref. [l7|, [HI for the case r]^ — {rj^). 

The KATRIN experiment, dedicated to measurement of electron neutrino mass from endpoint electron energy 
spectrum of /3-decay of tritium, can, in principle, also search for cosmic neutrinos via neutrino capture reaction p4l) . 
The KATRIN experiment (in construction phase) aims to measure rrii, with sensitivity of 0.2 eV [2^ using about 
50 fj,g of tritium corresponding to 5 x 10^* T2 molecules For this amount of target nuclei we find from Eq. ((39)) 
the number of neutrino capture events 



N^^ptiKATRIN) « 4.2 x 10-« ^ y-\ (40) 



Considering this estimate we conclude that observation of relic neutrino capture in KATRIN experiment looks rather 
unrealistic. Further comments on this subject will be given in sec. IIVI 

B. Neutrino capture on rhenium 

The cosmic neutrino capture by rhenium 

i^e ^e((5/2)+) ^ i«^Os((l/2)-) + e- (41) 

we analyze in the way similar to the previous case of tritium. The T-matrix determining this process coincide with the 
rhenium /3-decay matrix element shown in Eq. ()25|) . After neglecting the part of the electron wave function associated 
with emission of the Si/2-electron, which is very small, the squared and summed over polarizations T-matrix takes 
the form 

E \{f\T^'W = Gl BcTC'Re) F,iZ + 1,E,). (42) 
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Inserting this expression in Eq. (1331) . for the total capture rate per cosmic neutrino we obtain 

T^C^'Re) = ^-Gl Bore'' Re) F,{Z + l,E,) (43) 



with Ee = E., + E^- Ef. 

As in the case of tritium we replace 1/V with the neutrino number density rji, and get the capture rate per atom 
of rhenium 

r-(i87i?e) ^ }_Gl F,iZ + l,E,) Bare'' Re) Pe ^ M. (44) 

Substituting to this equation numerical values of nuclear matrix element from (|3ip and other constants we obtain for 
the capture rate per atom of ^^"^Re 

T'^C^^Re) = 2.75 x IQ-'''' y-\ (45) 

For a detector with 100 g of ^^"^Re, i.e. with 3.2 x 10^'^ rhenium atoms, we find the number of cosmic neutrino 
capture events 

N-^^, c 8.9 X 10-'^ y-'. (46) 

The MARE project will investigate the /3-decay of ^*^i?e with absorbers of metallic rhenium or AgReOi. It foresees 
a 760 grams bolometer. For this amount of rhenium the number of neutrino capture events is 



N^,p,{MARE) ~ 6.7 x 10-« ^ y-\ (47) 



The MARE detector technology can, in principle, be scaled up. With about 4 orders more rhenium the capture rate 
would become about two orders of magnitude larger than in case of the KATRIN experiment. Note that the KATRIN 
experiment can hardly be significantly scaled up due to safety limitations on the amount of highly radioactive tritium. 



IV. DISCUSSIONS AND CONCLUSIONS 



As was shown in the previous sections both with "^H and ^®^Re target nuclei the relic neutrino capture rate Eqs. 
(|40p - (l47p event numbers are extremely small and unobservable in the present and near future experiments. The 
gravitational clustering of relic neutrinos at the level of r]^/ {-q^) ~ 10'^ — 10* 0,[2l| can hardly change this conclusion. 
On the other hand this effect could be significantly stronger leading to a clustering about ?7iy/(?7iy) — 10^^ discussed 
in Ref. [s^. Then the capture rate becomes large reaching to about iV^^p^ ^ 10*" which looks very promising. 
However, a key point is distinguishing the relic neutrino capture signal from the /?-decay background. In this respect 
the main quantity to be analyzed is the ratio X^/Xp of the the partial rates of /3-decayA^ and relic neutrino capture 
Xu with the electrons in the energy interval around the endpoint Qp oi the /3-decay electron energy spectrum. Despite 
the /3-decay spectrum ends at Qp — and the electrons from the neutrino capture locates at + mu in a realistic 
experiment they overlap due to finite energy resolution A of a detector. As it was shown in Refs. [13: 13 neutrino 
capture events can be discriminated from the /?-decay background if m^/ IS. 2 or smaller. This condition is very 
challenging for experiments and is nearly independent of target nucleus and very weakly depend on the clustering 
rju/ iVu}- Only a strong clustering rj^/ {rji,) ~ lO^'^, considered in Ref. [13, may relax this condition and simultaneously 
significantly increase the total event rate of the relic neutrino capture. We note also that the considered searches for 
relic neutrinos via neutrino capture on /3-decaying nuclei even been unable to detect this process may set useful limits 
on the actual relic neutrino density iiv / {r]u) in the Earth vicinity. These possibility requires additional study. However 
it should be pointed out that there is a big uncertainty in this sort of studies due to unknown absolute mass scale of 
neutrino. If the forthcoming /3-decay experiments @,I23|, [25| set the value of m^^ within their sensitivities then the 
energy resolution necessary for successful searches of relic neutrinos will be within reasonable values like A ^ O.leV 
[iTj . which could be achieved in future experiments. 

In conclusion, we carried out a detailed analysis of nuclear physics and kinematical aspects of single nuclear /3- 
decay and relic neutrino capture on /3-radioactive nuclei. We focussed on '^H and ^^''Re isotopes to be used in the 
forthcoming KATRIN and in the planned MARE experiments, respectively. 

As to the prospects of relic neutrino direct detection through neutrino capture on /3-decaying nuclei we concluded 
that they are rather pessimistic with the present experimental methods and techniques. Even taking into account the 
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gravitational clustering, which greatly enhances the relic neutrino number density in galaxy clusters, does not improve 
this situation: the detection rate remains very small and hardly observable in the near future. However we noted 
that scaling MARE experiment up to several hundreds of kilograms of rhenium would offer an event rate significantly 
larger than in KATRIN experiment. This possibility, however, remains technically very questionable. 

Despite our rather pessimistic conclusions we point out that searches for cosmic neutrino capture via beta decaying 
nuclei deserves to be carried out. After all this is the only known direct way to probe this component of the universe 
and from non-observation of the capture process to extract limits on the local neutrino number density rj^, which could 
be complementary to the known cosmological and astrophysical limits. On the other hand one can never exclude that 
we are passing through a neutrino clump with a density much higher than predicted by the gravitational clustering 
[iM HH; • This would eventually make observation of cosmic neutrino capture realistic in even the forthcoming 
experiments. 
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Appendix: Partial electron wave functions 



Here, for convenience, we show several terms of the partial wave decomposition of the relativistic electron wave 
function in the Coulomb field of a uniform charge distribution of a nucleus: 

^elr, E) ^ ^-^1/2 (r, E) + (r, E) + ^p^/^ {r,E) + ... (A.l) 

In our notation we follow Ref. [27J. Keeping the leading terms in r we have for the partial wave functions: 
Si/2 wave: 

rs'^'ir,E)=( , \' 



Pl/2 



wave: 



P3/2 wave: 



VFUZ^E) u,iP). (A.2) 



-/-I f^-r Xs 

i ^ VFoiZ,E) jol-r UsiP). (A.3) 



g-2 [3f ■ p- (T r (7 ■ p] Xs 
f+2 [3f pcr p-(T r] Xs 

VFiiZ, E) fr-p+i7-r 7-pj Us{P). (A.4) 



The relativistic Fermi function Fq{Z, E^) {Fi{Z, Ee)) takes into account the Coulomb interaction of emitted Si/2 and 
Pi/2 {P3/2) electrons with the nucleus and is given, for instance, in Ref. (27| . The Dirac spinors are normalized as 
wt(P)u(p) = 1. 
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